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The spectral density of random matrices is studied through a quaternionic generahsation 
of the Green's function, which precisely describes the mean spectral density of a given matrix 
' under a particular type of random perturbation. Exact and universal expressions are found 

in the high-dimension limit for the quaternionic Green's functions of random matrices with 
O ■ independent entries when summed or multiplied with deterministic matrices. From these, 
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the limiting spectral density can be accurately predicted. 



Ph. I. INTRODUCTION 



A central problem in random matrix theory (RMT) is, and always has been, to determine the 
distribution of eigenvalues of a random matrix ensemble. For a given N xN matrix X the statistics 
• of the eigenvalues {A^-"'^^}^^ are captured in the spectral density 



' a generalised function of the complex variable A. Of particular interest are those ensembles for 

which the spectral density converges to a non-random limit as — )• oo, even more so when this 



limit is universal in the sense that it is independent of the distributions of the entries of the 
matrix. For example, in the case that matrices {^at} are taken to be Hermitian with otherwise 
independent entries satisfying a Lindeberg-type condition, the so-called 'simple approach' of Pastur 
and collaborators provides a robust and straightforward framework to prove the convergence 
of q{X; Ai\f+Di\i) for suitable deterministic Hermitian matrices {Di\i}. The limiting spectral density 
in this case is given by the celebrated Pastur equation [5|. 

Unfortunately, it has long been recognised that many of the techniques used in the analysis of 
Hermitian random matrices simply do not apply to their non-Hermitian counterparts (see, for 
example, comments in loj] and [s]), and universal results like the Pastur equation are notoriously 
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hard to come by. A case in point is the famous Circular law which, after many years of research 



and several partial results 



6l4lll|. has only recently been proved in full generality 



The purpose of this paper is to show how the methods of the simple approach can be re-engineered 
for use in the general setting of non-Hermitian matrices. In doing so, we will prove an extension of 
the Pastur equation with no Hermiticity requirements placed on either the random or deterministic 
part. A similar result is found for the product of a random and a deterministic matrix, a result 
which (for obvious reasons) has no analogue in Hermitian RMT. 

Amongst the numerous physical applications of non-Hermitian random matrix theory, we mention 
only a pedagogical example of quantum chromodynamics (QCD) with finite chemical potential. 
The density of states, derived first by Stephanov in [isl], can be recovered easily as an example of 
the sum rule. 

In common with many techniques in the study of Hermitian random matrices, the simple approach 
relies upon the well-behaved nature of the Green's function of Hermitian matrices, defined by 



N 

t 

For A away from the real axis (to which the eigenvalues of Hermitian matrices are confined) the 
Green's function is analytic and its imaginary part gives a smooth and A^-independent regularisa- 
tion of the spectral density. If X is non-Hermitian however, the eigenvalues invade the complex 
plane and the Green's function provides no such regularisation. Instead, we have the exact relation 

q{\-X) = --^G{\-X), (1) 

where d/d\ is the anti-holomorphic derivative. It has been frequently suggested to study the 
spectral density of non-Hermitian matrices through a regularised form of the Green's function or 
some related object^. This is usually achieved by association with an Hermitian proxy, of which 
there are several equally good (and often equivalent) choices. In an influential series of papers 



15l4l7l| Feinberg and Zee worked with 2A^ x 2A^ block matrices of the form 



{X - A)t ie 
a process they christened 'Hermitianisation 



H=\ / ' h (2) 



^ For instance, an electrostatic potential introduced in [l^ and used frequently thereafter. 
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Around the same time, Janik, Nowak and collaborators proposed a similar block extension tech- 
nique, obtaining a generalisation of the Green's function with a quaternionic structure [isl. [l^. In 



these and subsequent works 



2C 



the application of free probability theory to this quaternionic 
formalism has yeilded many interesting results, including for sums of unitary random matrices [2^ 
and infinite products of large random matrices 23] ■ In this paper, we will take a different direction 



and instead seek to adapt the techniques of the simple approach to Hermitian RMT, obtaining 
results that hold for matrices with independent entries of unspecified distributions. 

Before we procede, it should be noted that there is a significant drawback to the premise of 
introducing an e > regularisation to the Green's function of a non-Hermitian matrix. Simply 
put, if the matrices involved are not normal, there may be parts of the complex plane far from the 
spectrum in which the Green's function is nevertheless very large. In practical terms, this causes 
great difficultly in justifying the exchange of the limits N ^ 00 and e — )• 2J|. One route around 



this problem involves the analysis of least singular va 
a large part of recent work on the Circular law j^, 0, 1^, 1^ 



ues of the matrices involved, and has formed 



251 ]. A tentative link is often made to 



the influence of the pseudospectrum, though this idea is rarely expanded upon. 

In the present paper no attempt is made to explicitly tackle this problem, however, we are able to 
offer a remarkable relation between the (regularised) quaternionic Green's function and the mean 
spectral density of a given matrix under a particular type of random perturbation. 

The main results of the paper are stated in Section 2, together with some brief examples of the 
sum and product rules and a conjecture regarding a certain class of random matrices with an 
asymptotically spherical spectral density. Proofs of all the main theorems are provided in Section 3 
and the final section contains a discussion of the limitations of the work and some possible directions 
for future research. The notation used to define and manipulate the quaternionic Green's function 
is introduced in the remainder of this section. 



4 



Notation 



We will be doubling the size of our matrices. To simplify the formulas, we introduce the following 
notation: 

Denote by X the 2A'' x 2A'^ matrix composed of A^^ blocks of size 2x2 whose block is given 



We will always use Roman indices i, j, k, I to refer to the 2x2 blocks of boldface matrices, rather than 
the individual entries. All sums over Roman indices run from 1 to AT. To work with quaternions, 
we introduce the 2x2 matrix 



Now, if g = a + 6j is a quaternion (i.e. a and b are complex numbers and j is a quaternionic basis 
element), then we have the matrix representation 



This is an isomorphism, and |g| = ||q||, where || • || is the spectral norm. When q is a 2 x 2 matrix, 
and X a 2N x 2N matrix, we use the shorthands 

qX = [q® In)X , and {X + q) = X + q ^ In ■ 

In addition to the usual operations, we define an elementwise product for quaternions 

(a + bj) ■ (c + dj) = ac + bdj . 

Note that the matrix representation of an elementwise product of quaternions is not the same as 
the usual elementwise product of the matrices, in fact we use 



by 






(3) 
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II. MAIN RESULTS 



Let X be an X matrix, A a complex variable and e a strictly positive real number. Putting 
q = X + ej , we define the 2N x 2N 'resolvent' 

n{q;X) = {X-qr' . 

To connect with other approaches, note that there exists a permutation matrix P such that 

{X-q)=p\ HP-' , 

where H is the 'Hermitianised' block matrix given in ([2]). The quaternionic Green's function of X 
at q is then defined to be the quaternion Q{q;X) with matrix representation 



N 

t 

Without the hypercomplex part, the quaternionic Green's function agrees with the usual Green's 
function, 

g{X + Oj;X) = G{X;X) + Oj. 

Adding a positive real regulariser e > 0, we apply ([1]) to obtain from Q a regularisation of the 
spectral density, 

Qs{X;X) = --Re-^g{X + ej;X). 
TT dX 

Interestingly, this regularisation is precisely the mean spectral density of X under a particular type 
of random perturbation. 

Theorem 1 (Perturbation Formula). Let X be an arbitrary N x N matrix and e a strictly positive 
real number. Suppose A and B are random N x N matrices, with independent standard complex 
Gaussian entries, then 

E g{X] X + eAP-') = ge{X; X) . 

A short proof is presented in the next section, using a matrix generalisation of the Mobius trans- 
formation. As a corollary of Theorem 1 we see that mean spectral density of the matrices AB~' 
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is, regardless of their size, the uniform distribution on the Riemann Sphere^. We are immediately 
prompted to ask if this result holds in the limit N ^ oo for any distribution of the entries of 
and Bn, in an analogue of the circular law: 

Conjecture (The Spherical Law). Let {^at} and {i?Ar} be sequences of matrices of independent 
complex random variables of zero mean and unit variance. Then the spectral densities of the 
matrices A^B^^ converge to the uniform density on the Riemann sphere. 

This phenomenon was also noticed by Forrester and Mays [2^ and can be easily derived in a non- 
rigorous fashion using the techniques developed below, however a full proof is likely to require a 
more in-depth analysis. 

The other main results of this paper concern the quaternionic Green's function itself, and specifi- 
cally the ease with which universal predictions can be made about the limiting regularised spectral 
density of sums and products of random and deterministic matrices. 

Suppose we are in possession of an infinite array of complex random variables with joint 

probability space (ri, J^, P). We assume the ^jj to have the following properties: 

Al) K^ij = for all i,j 

A2) El^ijf = 1 for ah i,j 

A3) There exists a finite constant such that E < for all i,j 

A4) All S^ij are independent, except for the covariance E^jj^j, = r, where r G [0, 1]. 



A normalised N x N random matrix can then be constructed by taking 

1 



(4) 



Introduced in 



14}], the parameter r controls the degree of Hermiticity of vliv- At r = 1, we have that 
Afyf is an Hermitian Wigner-class matrix, whilst at r = the entries are completely independent. 
In our calculations, r will only appear as the real part of the quaternion t = t + j. 



In fact, the full jdpf of eigenvalues for matrices of this type was found recently in [26l |. 
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In the results below, we will characterise the quaternionic Green's functions of the sum or product 
of such matrices with deterministic matrices {-D^v} satisfying some or all of the assumptions 

Dl) The quaternionic Green's functions of {-Cat} converge pointwise to the limit Qd 
D2) The quaternionic Green's functions of {D^^} converge pointwise to the limit Gd-^ 
D3) There exists a constant d G [0, oo) such that sup ||^Ar|| < d and sup < d. 

The last point is a technical assumption made for the sake of simplicity and, as we note in a later 
example, may not be strictly necessary. 

Theorem 2 (Sum Rule). Let {A^} be a sequence of random matrices given by and let {D]\f} 
be a sequence of deterministic matrices satisfying Dl. Fix a quaternion q = X + ej, where e > 1. 
Then 

g{q;DN + AN)^Q{q), 

where Q[q) satisfies 

Q{q)=QD(q + t-Q{q)), (5) 

Theorem 2 is a straightforward generalisation of the Pastur equation for the sum of deterministic 
and random Hermitian matrices, indeed at r = 1 and e — t- 0, equation ^ precisely is the Pastur 
equation. For the case r = 0, an equivalent result has already been found using potential theory 

m- 

Theorem 3 (Product Rule). Let {^at} be a sequence of random matrices given by ^ and let 
{-Djv} be a sequence of deterministic matrices satisfying D1-D3. Fix a quaternion q = X + ej, 
where e > 2d. Then 

g{q; DnAn) a g{q) = -{t- qY^Qb ( - g (t • a)"') , (6) 

where Q satisfies 

§=-q-'GD-^{-{t-§)q-'). (7) 

Unlike Theorem 2, this result is not related to any in Hermitian RMT for the simple reason that 
the space of Hermitian matrices is not closed under multiplication. 
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Examples 



The statement of the sum and product rules given above concerns the behaviour of the quaternionic 
Green's function in the hmit N — )• oo, for a fixed regulariser e, which is taken to be large. In light 
of Theorem 1 we are in effect computing the limiting spectral density of matrices under a large 
perturbation. However, as the following examples will demonstrate, accurate predictions about the 
spectral densities of sums and products of matrices satisfying conditions A1-A4 can be found by 
naively taking £ = in the final equations. 

The Elliptic Law: 

The well-known elliptic law occurs naturally. Taking either = in the sum rule, ov = In 

p 

in the product rule, gives Q{q;Aj\[) — > Q{q), where 

-1 



Q{q) = -[q + t-Q{q) 

Writing Q{q) = a + we send e — )• and assume that /3 stays strictly positive, obtaining 

(a + /3i)(A + ra + /3j) + l = 0. (8) 

The support of the spectral density is restricted to the region allowing a solution with /3 > 0. The 
hypercomplex part of ([8]) gives 

( X y 
a = — \ h i- 



T + 1 T-1 

where \ = x + iy, and the complex part gives 



/5 = A 1 



> 2 / \ 2 

X \ y 



The condition /3 > determines the elliptic support, and taking an anti-holomorphic derivative 
yields the spectral density inside that region: 



pW = { 



^ when (^)'+(^)'<l 



otherwise. 

V 

Sum Rule: 

Let us consider the following random matrix model for the Dirac operator in QCD with finite 
chemical potential, considered in |13;], 

ztItv + M . 
D=\ . n , (9) 

iA\ + fi 
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where 747v is drawn from the Gaussian Unitary Ensemble and fi > is the chemical potential. 
Multiplying by an imaginary unit and filling the diagonal blocks, we suppose that in the limit 
— )• oo, the spectral density of (jO]) may be recovered from that oi A]\f + iM]\f, where 

Mn = 

For the quaternionic Green's function, the sum rule states that Q{q;AN + iM^) a + f5j, where 
equation ([5]) now reads 

a + l3j = ]^{iix-q-a-l3j^ + q + a + /3j^ . (10) 

Assuming the limit e — )• 0, we write q = X. The support of spectral density is then given by the 
region allowing a solution of (jlOp with /3 > 0, in this case determined by the condition 



l-lx'-\y'[l + 2{y' - /.2)) ' (y2 - f,') + {y' - /.^) > . 

Inside this region, one may solve for a and take the anti-holomorphic derivative to determine the 
spectral density 

1 / y2 + ^2 



n 

As expected, this is precisely the density first recovered by Stephanov in [13;], rotated by 7r/2. In 
that work, only the mean of the spectral density was computed, however the sum rule tells us that 
the quaternionic Green's function converges in probability as — )• oo, suggesting weak convergence 
in spectral density. Moreover, only Gaussian distributed random matrices were considered in 
[3], whereas the sum rule predicts the limiting density to be universal in the sense that it is 
independent of distribution of the entries of ^at. An analogue of the construction ([9]) for orthogonal 
and symplectic ensembles has been studied numerically in [3], where, interestingly, the limiting 
densities were found to be different to those for the unitary ensemble. 

Product Rule: 

We compute the limiting spectral density for the product matrix D^An, where the Aj^i are given 
by dlj) with T = 0, and Dj^j is a diagonal matrix with entries Da drawn independently from the 
standard Cauchy distribution^. With this choice oi D^, we have the limits 



1 r°° 1 

g(q-Dr,)^gDiq) = - / -—^ir-q)-'dr., 
vr 1 + H 



^ Notice that although this choice fails the technical assumption D3, the result still appears to hold. 
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FIG. 1: Blue histogram: the log-moduH of eigenvalues of a single random matrix DA of size N — 10,000, 
where the entries of A are independent complex Gaussians with variance 1/A^, and Z? is a diagonal matrix 
of Cauchy random variables. Black line: the density predicted by the product rule (Theorem 3), given in 
this case by equation ([T2|) in the text. 



and Q{q;D^^) — )• Gd{q), as N ^ oo. As before, we assume the e — )• limit, taking ^ = A. Then 
Q{q) = a + where the product rule ([7]) gives 



vr ./-no 1 + 



+ /3j dr . 



Performing the integral and solving for /?, we obtain /3 = ( — |A| + + 4) /2. Returning to ([6]), 

we reach 



N 



^ /"OO 

TT 



1 



1 + r2 



-1 

A + r/3j ) dr 



2A 



|A|2 + |A|vW + 4 



and finally an expression for the limiting spectral density. 



1 



p{X) = - 



^ V|Ap + |A|V|A|2 + 4 |A|2 + |A|V|A|2 + 4 + 4y ' 

To provide an effective comparison with numerical data, we change variables to 7 = In |A|, whose 
distribution is given by the pdf 

e'fVe'^-y + 4 - e^T 



^(7) = 2- 



(12) 



e^Ve^T + 4 + 6^7 + 4 ' 

Figure [J shows a histogram of the log- moduli of the eigenvalues of a single such random matrix of 
size N = 10,000 alongside the predicted density 2^(7). 
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III. PROOF OF MAIN RESULTS 



A. Proof of perturbation formula 



Proof of Theorem 1. For the case N = 1, the ratio of two standard complex Gaussian random 
variables takes the uniform density on the Riemann sphere. One proof of this fact comes from the 

observation that the density generated is invariant under a class of Mobius transformations. We 
generalise this idea to matrices. Begin by noting that, for an arbitrary N x N matrix X, 

ee{X;X) = -lTrA((x - A)t(X - A) + e'y\x - A)t , 

and 

g(X;X + eAB-^) = -—Ti^ieAB-^ + X - X)~^ . 
Theorem 1 will then follow from the stronger claim that for any matrix Y 

A little rearrangement leads to the equivalent statement E/(^S~^) = 0, where we have introduced 
the matrix Mobius transformation /, given by 

Notice that if | AB'^, then F 1^ j iH- /(AS"^), where F is the IN x IN block matrix 

Vt -e 
e Y 



It will be useful to 'normalise' /, introducing 

f{Z) = (yty + e") {Y^Z - e) {eZ + Y) (yyt + e^) V2 
Then the normalised form of F is given by G SU(2A^), 



F 



/ \ -1/2 / 

y YY^ + e^ j V ^ ^ 
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Now, let A and B be independent complex Gaussian matrices, with joint pdf 



P 



A 
B 



TT 



exp < — Tr 



At St 



A 
B 



and let 
B' is given by 



Ia'\ 












\b'I 







. We perform a change of variables to find that the joint pdf of A' and 



2Af2 



TT 



exp < — Tr 



A^ St ) (Fti?) 



,-1 



A 
B 



detF 



-1 



A 
B 



Thus the density p is invariant under multiplication with F and we can conclude that the distri- 
bution of the random variable AB~^ is invariant under /. 

So to prove that Ef{AB~'^) = 0, it will suffice to show that EAB^'^ = 0. But A and B are 
independent, = 0, and it can be shown (see Q) that = VNn < oo, so we are 

done. □ 



B. Preliminaries for working with the quaternionic Green's function 

For the proofs of Theorems 2 and 3, a number of standard tools will be of repeated use. 

The resolvent identity. For N x N matrices A and B, and quaternion q, we have 

n{q; A) - n{q; B) = n{q; A){B - A)n{q; B) . (13) 
This is a consequence of the more general expression for any X and Y , 

X-^ - Y-^ = X-^{Y - X)Y-^ . 

The resolvent bound. For any matrix A and quaternion q = X+ej, we have the following bound 
on the norm of the resolvent and its blocks 

\\n{q;A)ij\\<\\n{q;A)\\<y (14) 

Again, this is a special case of a more general result - for any matrices X and Y, with Y 
invertible, we have 

1 



|(X + rj)-.i|| < \\{X + Yj)-^\\ < 



13 



where x and y are the smahest singular values of X and Y . 



The Cauchy-Schwarz inequality. Let {xi}f^^ and {yi}f^i be collections of 2 x 2 matrices. We 
have the following analogue of the Cauchy-Schwarz inequality: 

2 



E 



< 



(15) 



Integration by parts. Let / : C" — t- C be a continuous and differentiable (but not generally 
analytic) function, with bounded second order partial derivatives 



sup 



d 



dxidxj 

Then if F = f{Cu, • • • , Cnn), we have 



d 



dxidxj 



fix] 



d 



dxidxj 



<Cf. 



dF 



dF 



dF _ ^ dF 



(16) 



where |A;| is bounded by some constant depending on and Cj. Proof is by Taylor's 
Theorem. 



In addition to these general facts, the main part of the work in proving Theorems 2 and 3 comes 
down to the application of two central results. 

Lemma 1. Let An be a random N x N matrix given by and let X , Y and Z be arbitrary 
deterministic matrices of the same size, with Z invertible. Define 

R={XAn + Y + Z3)-^ , and G=^Y.^i- 

i 

Then 



E 



G-EG 



< CN-^ , 



where C is a constant depending on X,Y and Z. 



This is simply a slightly more general version of the statement that the quaternionic Green's 
functions of the matrices we are interested in are self- averaging in the limit N ^ oo. This property 
is crucial if we are to extract useful information about the behaviour in that limit. The proof is 
based on an elegant martingale technique from a paper on spin-glasses 3l|. 
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Proof. Let P = : 1 < i < j < N}. We label these pairs by introducing the bijective 

numbering -H- p where p G {1, \P\}- The names and p will be used interchangeably; 
the meaning should be clear from the context. Let = T and recursively define the sub-a-algebras 
•^P — o'l-T^-i, '^jj, Cjj}- Introduce the martingale 

Ap = E (G I Jp) - E (G I , 



so that 

\P\ 

= G-EG. 

p=i 

We plan to bound each Ap, to do so, we consider the fictitious situation in which the blocks (A^r) 
and (Atv) •• are removed. Write A^^'' for the matrix obtained from hy setting = = 0. 



Introduce 



E^^^^ = {XA^KY + Zjy\ and G^^^') = 1 J] . 

k 



The resolvent identity (fT3]l provides 

R = R^'i^ - RX [An - A%'^) 

and thus 



where the error term Kij is given by 



kk 



The Cauchy-Schwarz inequality (jl5p provides a bound for -K^ij, since, for example 

2 

">:(i?x),^(i?x)t 



A; 



< 



^jk ^Jk 



< 



CM, 



N 
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where Mij = max{|^jj|, and C is the constant coming from the resolvent bound (jl4p . We 

can conclude \\Kij\\ < CkM^N"'^/'^ , for some constant Ck-, and thus 

|A„|| = E (G I Jp) - E (G I Jp_i) 
E {Eij I Jp) - E {Eij I 
<Ci^7V-3/2(M,,--EMi,). 
Burkholder's inequality then gives, for some constant Ca, 



E 



G-EG 







3 




E 




< CaE 






p=i 







3/2 



3/2 



where Ca/ < 2Cg + 2 is a bound for EM?-, and we have repeatedly used Jensen's inequality. The 
desired bound is then given by 

2 



E 



G-EG 



< CN-^ , 



where C = {ACaC^^Cm?^^ ■ 



□ 



With self- aver aging established, we next require a mechanism by which we can convert the general 
statement of the resolvent identity ()13p into an equation for the mean of the quaternionic Green's 
function. 

Lemma 2. Fix a quaternion q = X + ej, with e > 1. Let A]\f be a random N x N matrix given by 
^ and let X be an arbitrary deterministic matrix of the same size. Define 



R= {XAn - q) ^ and G = N'^"^ (RX) .. 



Then 



E (Ajv-R) = - (E (t • G)) (E R) + Kn 
where t = t + j, and \\Kiq\\ as N ^ oo. 
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Proof. We compute a generic block 




(17) 



Notice that the entries of R depend continuously upon the ^, and the resolvent bound gives a 
bound on the first and higher order derivatives. We are thus able to apply the integration by parts 
formula ([16]). The derivatives are given by 



dR 1 „^ / \ „ ^ OR 

d^al3 VN \ 1 / 



1 / 1 

RX\ I l^^R, 

N loo 



where lai3 is the 2N x 2N block matrix containing a copy of I2 in block (a, /3) and zeros elsewhere. 
Applying this to (fT7|l . we thus find, after some tedious algebra, 



j 

= -E ((t • G)R,k ) - ^ E (^^^- • (^^),.) ^^•'^ ' 



where Sij = ^^fj + (E^ij^jj)j. Notice that assumptions A1-A4 imply a universal bound for \sij\, 
and the Cauchy-Schwarz inequality together with the resolvent bound give a constant C such that 



1 

N 



Finally, we split the expectation 



< CN- 



E{{t-G)R) - (E(^•G))(Ei^) 



as N ^ 00, 



since E IIG — EGII — >• by Lemma 1, \t\ < 00 and ||-R|| is bounded. 



□ 



C. Proof of sum and product rules 



Proof of Theorem 2. Fix g = A + ej with e > 1. We write the shorthands R^ = 72.(g; + ^at). 
Gat = Q{q;DN + An) and Qn{<i) = Q{<1]Dn + ^at)- Now, applying the resolvent identity and 
Lemma 2, we obtain 

Ei?Ar - (£>7v - q)-^ = -{Dn - Q^^EAnRn 

= (Dn - q)-' (E (t • Gn)) (E Rn) + Kn , 
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where H-K'ivll — )• as — )• oo. Rearranging, we have 

ERn = (^Dn -q-t-EGN) \k'n, 

where the resolvent bound gives — >• as — )• oo also. Summing over the diagonal blocks, 
we deduce 

EgNiq)=g{q + t-EgNiq);DN)+kN, (18) 
with |A;Ar|— )-OasA'"— T-oo. Define the functions 

fN{g) = Gdn (q + t- g) +kN , f{g) = QD{q + t- ■ 

Since e > 1, the resolvent bound gives that each /^r is a contraction with parameter e~^, and thus 
the pointwise limit / is also; it is therefore continuous, with a unique fixed point which we call 
Q{q). Finally, Lemma 1 gives E \QN{q) - QN{q)\'^ ^ as iV ^ oo, and we conclude from (llSp and 
Tchebychev's inequality that 

GN{q) -^G{q)- 

□ 



Proof of Theorem 3. The proof is very similar to that of Theorem 2, so we give only the main steps 
of the derivation. 

Fix q = \ + ej with e > d. Again we write the shorthands Rn = 7t{q; DnAn) and GN{q) = 
Q{q; D]\fAj\f). Let Gn = N^^ J2i{^NDj\f)ii. Applying the resolvent identity and Lemma 2, we 
obtain 

ERn + = q^^DnEAnRn 

= -q-^DN{E it-GN))iERN)+KN, 

Rearranging, we have 

ERN = -{E{t- Gn)) [Dn + q{E{t- Gn)) '^y\ K'^ , 

where the resolvent bound gives H-f^yvll — as A^ — )• oo also (it is here that assumption D3 and 
the requirement e > 2d are used). Summing over the diagonal blocks, we deduce 

EgNiq) = -{E{t-gN)y^g(^-q{E{t-gN)y^;DN) +kN, (19) 
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where Qn '^s the quaternion with matrix representation Gn, and \kj\f\ — )-OasA^— t-oo. The 
equation for G is found similarly; let Rn = RnDn, then 

ERn + q'^DN = q-^Dn^AnRnDn 

= -q-^Dn (E (t • Gn)) (E Rn) + Kn ■ 

and thus 



E 



Gn = -q-'G [-{t ■EgN))q-';D]^'^ + kN , (20) 



where as usual |fcAr| — )-OasA^— )-oo. The proof is completed in the same fashion as Theorem 2, 
with equations (fT9|) and (pOj) providing ^ and ([7]), respectively. □ 



IV. DISCUSSION 



The purpose of this pap er has been to marry the simple approach to Hermitian RMT [Jll to 



the tricks used in 15142 ll| to deal with non- Hermitian matrices, thereby obtaining techniques with 
which to handle sums and products of random and deterministic matrices. As shown in Theorem 1 
the resulting theory in fact applies to the mean spectral density of such matrices under a particular 
type of random perturbation. In practice it appears, as evidenced by the examples, that Theorems 
2 and 3 can indeed be used to predict limiting spectral densities in the absence of this perturbation, 
though this aspect of the theory has not been rigorously proven. There are two main difficulties 
in the justification of the exchange of limits N ^ oo and e — t- required to make the approach of 
the examples rigorous. 

First, in the statement of both results we assume a minimum size for e with the purpose of making 
unique fixed points of equations ([5]) and ([7]) easily available. The same trick is used in the analogous 
theory of Green's functions of Hermitian matrices, however, in that case it is easily justified by 
appealing to analyticity; one simply determines the Green's function far from the real axis and 
uses analytic continuation to return. Things are not quite so straightforward in the quaternionic 
case. A similar argument is still possible if one notes that for fixed X, g{X + £j;X) is a ratio of 
polynomials in e and that the zeros of the denominator are confined to the imaginary axis. One 
may feasibly then promote e to a complex variable on a strip containing the real line and apply 
analytic continuation as before. The only drawback here is conceptual; preserving the analyticity 
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of ^ in e effectively destroys the quaternionic analogy since we must work with and not |ep, as 
would result from the matrix representation ([3]). 

The other problem is more fundamental. For a sequence of matrices {X]\f}, the convergence 
in probability of G{X + £j;X]\f) translates to the weak convergence of g£{X;X]\f) and does not 
necessarily reveal anything about the limiting behaviour of the unregularised densities g{X]X]\[). 
If X]\f happens to be normal, then it is straightforward to compute 

Iff £ 



0s{X;Xn) = - {^r~\ ri2 0{i^;XN)dfi . 

^ Jc + Im - Ary 

Since, in this case, the smoothing of Qs is independent of N, the weak convergence of g{X;Xfyf) 
follows easily. Without normality the same cannot be said, in fact it is easy to construct examples 
for which lim£_j.o hniAf-i-oo Xj\[) and limAr^oo q{X; Xj\i) are entirely different*^. 



As mentioned in the introduction, this issue is by no means new or unique to the quaternionic 
Green's function. Some authors have treated the problem carefully 0,0,0, 10, 12, 25], usually by 
techniques that involve the bounding of the least singular values of the random matrices involved. 
Such methods may well be adapted to prove the convergence of densities q{X;Dn^ + An) and 
g{X; Dj\fAj\f) to the limits predicted by the sum and product rules, thereby completing the work of 
the present paper. 

We should point out however, that Theorem 1 suggests the regularised density, for small e or not, 
is an interesting and potentially useful object which is well worthy of study in its own regard. 
This result also offers a strong heuristic argument for the correctness of the techniques used in 
the examples - if we are dealing with very large and fully random matrices, the addition of an 
infinitesimal random perturbation should not change the spectral density. 

There is some room for improvement in Theorems 2 and 3, both in relaxing the conditions and 
strengthening the results. For the sake of simplicity and brevity, we assume in A3 the bound 
El^ijf < which could almost certainly be dropped in favour of a weaker condition, or possibly 
forgotten entirely as in 3]. In a similar vein assumption D3 may be extraneous, as suggested by 
the example. Lastly, it is possible that the convergence in probability of G may be traded up for 
almost sure convergence, which in turn would provide strong rather than weak convergence of g^, 
however this may well require entirely different methods. 



Banded Toeplitz matrices are a good choice. 
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The final opportunity for future research worth mentioning is the conjectured 'Spherical law', a 
proof of which, by any method, would be very interesting. 
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